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Spherically symmetric space-times provide many examples for interesting black hole solutions, 
which classically are all singular. Following a general program, space-like singularities in spherically 
symmetric quantum geometry, as well as other inhomogeneous models, are shown to be absent. 
Moreover, one sees how the classical reduction from infinitely many kinematical degrees of freedom 
to only one physical one, the mass, can arise, where aspects of quantum cosmology such as the 
problem of initial conditions play a role. 
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One of the main issues to be addressed by quantum 
gravity is the singularity problem of general relativity. 
While the classical theory is very successful in describ- 
ing space-time on scales which can be probed today, it is 
incomplete because it predicts the generic presence of sin- 
gularities: boundaries of space-time which can be reached 
by observers in a finite amount of time, but at which 
point the theory becomes inapplicable. Usually, curva- 
ture or energy densities and tidal forces diverge there, 
implying unphysical conditions. 

One explanation is that the picture of a smooth space- 
time underlying the classical theory is only appropriate 
at large scales, while at small scales the structure is dis- 
crete. This has indeed been substantiated by looking at 
cosmological models in loop quantum gravity 0, where 
the discrete quantum geometry has been shown to re- 
move singularities 0]. Even at the classical level there 
are indications for a break-down of the smooth picture 
at small scales: The BKL scenario provides a scheme 
for the general approach to a classical singularity by con- 
sidering dominant contributions to the field equations 
in this limit. It turns out that only terms with time- 
derivatives remain such that spatial points decouple and 
their geometries can be described by the most general 
homogeneous model. This Bianchi IX model is classi- 
cally chaotic 0] and so geometries in different points are 
completely unrelated, implying a complicated classical 
singularity with structure at arbitrarily small scales. 

Since the removal of singularities in loop quantum cos- 
mology applies for all homogeneous models in par- 
ticular the Bianchi IX model, one can combine this result 
with the BKL picture and expect that all singularities 
are removed by quantum geometry. However, the BKL 
picture has not been proven classically, and the above 
argument would require it to hold even in quantum grav- 
ity. In particular the latter point is questionable because, 
for one, already the approach to a single Bianchi IX sin- 
gularity in quantum cosmology is modified, removing 
the classical chaos 0, HI ■ It is thus necessary to study in- 
homogeneous models in loop quantum gravity and look 
at the singularity issue without assuming homogeneity. 
If this is possible, one can also test the validity of the 



BKL picture in the quantum context. This is what we 
will do here in the case of spherical symmetry, which is 
not only the simplest inhomogeneous situation but also 
allows interpretations for black holes. The same meth- 
ods apply to other models which do have local degrees of 
freedom, providing the first demonstration of the absence 
of singularities in inhomogeneous quantum gravity. 

Singularities. The main problem caused by a singu- 
larity is the fact that it presents a boundary to physical 
evolution. In order to see whether it persists in quantum 
gravity, then, the following steps have to be performed. 
This has to be done in a manner which is independent of 
coordinate or other gauge choices, and only potential sim- 
plifications resulting from the symmetry reduction should 
be used. One first has to locate classical singularities 
on the phase space of physical fields, the spatial metric 
q a b and extrinsic curvature related to q a b- Conditions to 
specify the singular part of phase space must be chosen 
such that any solution to the theory, which is a trajec- 
tory on phase space, intersects this singular part exactly 
when it develops a singularity. The solution space is in 
general quite complicated to study, but one can select a 
variable T on phase space which is transversal to the sin- 
gular part, a local internal time rather than coordinate 
time, i.e. which fulfills T = in a neighborhood around 
zero exactly at the singular part. Finally, one needs to 
write down the quantum evolution of geometry in the 
local internal time and check whether or not it stops at 
T = 0. If one can find a T such that the quantum evo- 
lution does not stop anywhere, the quantum system is 
non-singular. This is the analog of the classical notion of 
space-time completeness. 

We illustrate this scheme with isotropic cosmology 
where the phase space is 2-dimensional with the scale 
factor a (the spatial radius of the universe) and its time 
derivative. Singularities occur only if the scale factor 
vanishes such that a — specifies the singular part. An 
obvious local internal time (which in this case is global) 
is given by T — a, or with a slight modification the triad 
variable p with \p\ = a 2 and sgnp being the orientation 
of space. Using this variable makes no difference classi- 
cally, but is important in quantum geometry where triads 
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are basic variables. At the quantum level one can then 
first note that operators for p^ 1 are finite 0, indicating 
already that curvatures and energy densities do not di- 
verge, and most importantly that the quantum evolution 
is given by a difference equation for the wave function 
in p which does not stop at p = 2j. Thus, there is no 
singularity in isotropic loop quantum cosmology. 

Spherical symmetry. The case of interest here is 
spherical symmetry, where the kinematical phase space 
on which we have to locate singularities is infinite- 
dimensional and spanned by the metric components in 

ds 2 = q x (x)dx 2 + q Q (x)dSl 2 (1) 

(in polar coordinates) and their time derivatives. As an 
example we can look at the Schwarzschild solution for a 
black hole of mass M, q x = (1- 2M/x)~ 1 , q n = x 2 . The 
singularity is reached for x = 0, at which point both met- 
ric coefficients are zero. The question then arises which 
one, or both, of them must be zero as a condition for a 
singularity. It turns out that qa is zero only at the sin- 
gularity, while q x can also become zero elsewhere, i.e. at 
the horizon x = 2M, when one chooses a different gauge 
(e.g. with homogeneous coordinates in the interior). This 
point illustrates why gauge independence is essential in 
answering the singularity problem: even the very first 
step, finding where singularities would develop, depends 
on it. In fact, in this case we can choose our coordi- 
nates x and t at will, which affects the form of q x and 
points where it can be zero. In spherical symmetry, how- 
ever, the fact that qa = at the singularity is unaffected 
(even though, of course, gsi can change as a function of 
x when we change coordinates). 

We can now consider a spatial slice which locally, 
around a point xq, approaches the classical singularity 
such that qn(xo) — > 0. The above discussion shows that 
T = qn{xo) is a good local internal time, which com- 
pletes setting up the problem from the classical side. It 
now remains to formulate the quantum evolution in local 
internal time and check if it stops at T = 0. 

Quantum geometry. Again we first transform to triad 
variables \E X \ = qn and E v — ^/q x qn which become basic 
operators in quantum geometry. The (local) orientation 
of space around a point x$ is now given by $,giiE x (xo) 
where E x unlike E* can take both signs. Moreover, the 
discussion in metric variables shows that T = E x (xo) is 
our local internal time such that the situation, so far, is 
analogous to that in the isotropic case: triad variables 
lead to a local internal time which takes values at both 
sides of the classical singularity, T = defining a man- 
ifold in superspace rather than at the boundary. It is 
important to note that the introduction of triad vari- 
ables was seen as a necessary step toward a background 
independent quantization. Now it turns out that this 
also changes the singularity structure on phase space in 
a way which was important for removing cosmological 



singularities. Nevertheless, even though the singularity 
is now located in superspace, the classical evolution still 
stops there and is not able to connect from positive to 
negative T. This still has to be checked by the quantum 
evolution, the most crucial point. 

Quantum evolution follows from the Hamiltonian con- 
straint operator acting on states in the form of a lattice 
model with basis \k,j?) := ■ - « hn ff^" + \ V " ' 
where the integer labels fc e on edges are eigenvalues of 
the operator E x and the positive real labels fi(v) at ver- 
tices those of E v . Positions of vertices do not refer to a 
background space, and the lattice model represents the 
continuum theory. The constraint then acts by [Jjj 

H[N] k - . fc + - T, v N(v)(C (k) k - . fc + 
+C R+ (k) k - , k ++ 2 +C R -{k) k - . k +- 2 
+C L+ (k) k -+ 2 . fc + +C L - (fc) k — 2 . fc+ + • • ■ ) 

summing over all vertices of the lattice, the dots indicat- 
ing further terms such as a matter Hamiltonian whose 
detailed form is not important here. The known coef- 
ficients Ci(k) — Ci(k)Ci consist of functions C/(fc) of 
the edge labels and operators Cj acting only on the de- 
pendence on vertex labels pi. A general state is now a 
superposition \ip) — J2k p ^(^> A*)|k> P) whose coefficients 
i])(k,p) define the state in the triad representation. The 
constraint H[N]\ip) — has to hold true for all functions 
TV with independent values N(v), giving one equation for 
each vertex which in the triad representation takes the 
form 

Co (k)iP(k- ,k+) + C R+ (*)^(fe- ,k+ - 2) 
+C R -{k)Tp(k-,k+ + 2) + C L+ (k)i){k- -2,k+) 
+C , L _(fc)V(/c-+2,/c + ) + ■■• = () 

of a difference equation, where we have suppressed the 
vertex labels on which the Cj act and unchanged k. 

We now solve this set of equations with initial and 
boundary values for the wave function. To define a solu- 
tion scheme we proceed iteratively from vertex to vertex, 
starting at one side d of the lattice. We assume that the 
boundary values for all pig and k + (d) —: fc_ of the wave 
function as well as values for large positive k e — k$ and 
ko — 1 at all edges e are given, which means that we have 
specified the initial situation, e.g. by a semiclassical state 
specifying the initial slice far away from the singularity. 
The equation can then be solved for C R+ ijj(k-, k + — 2) 
in terms of values of the wave function specified by the 
initial conditions. This brings us one step further be- 
cause we now have information about the wave function 
at k + — 2 for a smaller edge label (our local internal time) 
evolving toward the classical singularity. 

Next, we have to know how to find tp from its image 
under Cr+. This can be done by specifying conditions 
for the wave function at small /x (which is not in the sin- 
gular part of minisuperspace but represents an ordinary 
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boundary) and happens in exactly the same way as in 
homogeneous models Before continuing, we notice 
that this indicates the presence of aspects of the BKL 
picture in quantum gravity. However, we still have to try 
to evolve through the classical singularity, i.e. k e = 0, 
which will be the main test. One crucial difference to 
cosmological models is that the coefficients C/(fc) are not 
only functions of the local internal time, k+, studied in 
the iteration but also of neighboring labels such as k- 
which do not take part in this difference equation but 
the dependence on which has been determined in iter- 
ation steps for previous vertices. This is clearly a new 
feature coming from the inhomogeneous context, and it 
has a bearing on the singularity issue. 

Singularities are removed if the difference equation de- 
termines the wave function everywhere on minisuper- 
space once initial and boundary conditions have been 
chosen away from classical singularities. The simplest 
realization is by a difference equation with non-zero co- 
efficients everywhere. However, this is not automatically 
the case with an equation coming from a general con- 
struction of the Hamiltonian constraint, and so has to be 
checked explicitly. Here, it turns out [lfj that a symmet- 
ric constraint indeed leads to non-zero functions C/(fc) 
which then will not pose a problem to the evolution. All 
values of the wave function, at positive as well as nega- 
tive k, are determined uniquely by the difference equa- 
tions and chosen initial and boundary values. The evo- 
lution thus continues through the classical singularity at 
zero k: there is no quantum singularity. Other quantiza- 
tion choices can lead to quantum singularities, providing 
selection criteria to formulate the quantum theory with 
implications also for the full framework. 

Consequences. We have shown that the same mecha- 
nism as in homogeneous models contributes to the re- 
moval of spherically symmetric classical singularities. 
Key features are that densitized triads as basic variables 
in quantum geometry provide us with a local internal 
time taking values at two sides of the classical singular- 
ity, combined with a quantum evolution that connects 
both sides. No new ingredients are necessary for inhomo- 
geneous singularities, only an application of the general 
scheme to the new and more complicated situation. 

As in cosmological models the argument applies only 
to space-like singularities such as the Schwarzschild one. 
The reason is that we evolve a spatial slice toward the 
classical singularity and test whether it will stop. A time- 
like or null singularity would require a different mecha- 
nism which is not known at present. Thus, cases like 
negative mass solutions seem to remain singular, which 
is a welcome property helping to rule out unwanted so- 
lutions leading to instability |l2j. 

This scenario and its form of difference equations does 
not only apply to vacuum black holes but also to spher- 
ically symmetric matter systems. In such a case, there 
would be new labels for matter fields, and a contribu- 



tion to the constraint from the matter Hamiltonian. As 
this does not change the structure of the difference equa- 
tion, the same conclusions apply. Moreover, models for 
Einstein-Rosen waves have a similar structure just with 
a new vertex label. Also in this case, with or without 
matter fields, the analysis goes through such that the 
absence of singularities can be demonstrated even in sit- 
uations with local gravitational degrees of freedom. 

There are differences between homogeneous models 
and these inhomogeneous cases, and the inhomogeneous 
analysis is much more non-trivial. In homogeneous mod- 
els there are several ambiguities in the constraint opera- 
tor, and several choices lead to non-singular evolution. In 
more complicated situations such as those studied here, 
not all options remain available. In particular, we had 
to use a symmetric ordering of the constraint in order 
to have non-vanishing coefficients of the difference equa- 
tion. In homogeneous models one can also work with a 
version whose coefficients vanish right at the singularity. 
The evolution then still continues since the value at the 
classical singularity simply decouples and does not play a 
role for the evolution. Instead, one can use the behavior 
to find dynamical initial conditions 01 • This is also pos- 
sible here for evolution in local internal time, but then 
the decoupled value at fc_ = is not determined and 
in general needed for the wave function at other values 
of k + . The inhomogeneous evolution would thus break 
down, and this choice of constraint is ruled out. 

There is a difference in the constraint operator we used 
compared to a common expression in the full theory 
|l4j . This issue is visible only in inhomogeneous mod- 
els, and consists in whether or not the constraint cre- 
ates new edges and vertices, or just changes labels of 
existing ones. We chose the second possibility, which 
has already been considered as a modification in the full 
theory |l5|. There, it can better explain the presence 
of correlations at an intuitive level, but makes checking 
anomaly-freedom more complicated. The main problem 
of an anomalous quantization would be that too many 
states could be removed when imposing the constraints, 
leaving not enough physical solutions. This issue can be 
checked here with the constraint we used. If there is no 
matter field present we expect just one physical degree 
of freedom, the Schwarzschild mass M. In our solution 
scheme we started with a boundary state ipg correspond- 
ing to this degree of freedom, and with this state being 
free it is already clear that we do not lose too many states. 
It is even possible to check whether or not the number 
of independent physical solutions is correct, i.e. not too 
large either. In the iteration we solve one difference equa- 
tion for tjj at each vertex, such that any freedom here 
would provide new quantum degrees of freedom. Since 
the difference equation for ip has the same form as that 
in homogeneous loop quantum cosmology, the number 
of quantum degrees of freedom is formally related to the 
initial value problem of quantum cosmology. A possible 
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physical meaning is to be checked in explicit examples. 

In the isotropic case there are indeed dynamical ini- 
tial conditions following from the dynamical law [T3L Hif 
which, if realized in our context, would imply that solu- 
tions for ip are unique and the mass is the only quantum 
degree of freedom. However, these conditions rely on the 
fact that leading coefficients of the difference equation 
can vanish, which we have ruled out for inhomogeneous 
models. Moreover, the uniqueness of a quantum cosmo- 
logical wave function depends on the pre-classicality con- 
dition of . Other mechanisms to select unique cosmo- 
logical solutions are thus needed, such as from observ- 
ables or the physical inner product This issue is 

quite complicated for difference equations in particular 
in anisotropic models 0|. Nonetheless, a simple count- 
ing of free variables supports the connection to initial 
conditions: The vacuum spherically symmetric case has 
difference equations in three independent variables, an 
edge label k and two neighboring vertex labels /i. Homo- 
geneous loop quantum cosmology gives rise to an equa- 
tion of similar structure and also three variables, so if 
we assume that there is a mechanism for a unique so- 
lution it will also apply to black holes of a given mass. 
Adding matter fields (or more gravitational freedom as 
in Einstein-Rosen) increases the number of independent 
variables to five in inhomogeneous models (two new ver- 
tex labels) as opposed to four in homogeneous matter 
models. The type of difference equations thus agrees in 
homogeneous and inhomogeneous models in vacuum, but 
not when local degrees of freedom are present. 

Thus, the structure of the Hamiltonian constraint 
equation from loop quantum gravity can potentially pro- 
vide explanations for issues as diverse as the singular- 
ity problem in cosmology and black hole physics, initial 
conditions in quantum cosmology, the semiclassical limit 
and issue of quantum degrees of freedom. We emphasize 
that many of these connections still have to be checked 
in generality. Still, such connections between seemingly 
unrelated issues in quantum gravity can be seen as sup- 
port for the internal consistency of the whole theory and, 
hopefully, provide guidance for future developments. 

We can finally come back to the approach to a classical 
singularity and the BKL picture. Our results here do not 
rely on an extension of the BKL picture to the quantum 
situation. First of all, the situation is conceptually differ- 
ent because evolution is now studied for a wave function 
in local internal time T, rather than the spatial metric 
in coordinate time. Nevertheless, at first sight a similar 
picture arises here from the quantum equation: as used 
in the previous arguments, the equations can be reduced 
to ordinary difference equations in T, where neighboring 
edges just contribute via an inhomogcncity of the differ- 



ence equation. The inhomogeneous situation, however, 
does play an important role right at the classical sin- 
gularity where some versions which would be allowed in 
homogeneous models are ruled out. 

Given that the techniques necessary for the quantum 
theory are similar to lattice models, it is easy to imple- 
ment them in numerical quantum gravity. This opens the 
door to numerical investigations of many problems that 
are still actively pursued in classical gravity 0, such as 
the approach to classical singularities and the issue of 
gravitational collapse and naked singularities. This re- 
quires studying horizons in addition to classical singular- 
ities, which can also be done at the quantum level [lflj . 
As we have seen, there are many non-trivial quantum ef- 
fects which play together in just the right way to ensure 
the absence of singularities, which has prospects for other 
effects in the physics of black holes [20j . 
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